Dark Matter from Minimal Flavor Violation 
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We consider theories of flavored dark matter, in which the dark matter particle is part of a 
multiplet transforming nontrivially under the flavor group of the Standard Model in a manner 
consistent with the principle of Minimal Flavor Violation (MFV) . MFV automatically leads to the 
stability of the lightest state for a large number of flavor multiplets. If neutral, this particle is an 
excellent dark matter candidate. Furthermore, MFV implies specific patterns of mass splittings 
among the flavors of dark matter and governs the structure of the couplings between dark matter 
and ordinary particles, leading to a rich and predictive cosmology and phenomenology. We present 
an illustrative phenomenological study of an effective theory of a flavor SU{3)q triplet, gauge singlet 
scalar. 



The observational evidence for dark matter (DM) pro- 
vides one of the few solid empirical motivations for new 
particle physics Indeed, the simplest explanation of 
DM is a new massive particle that is stable on cosmologi- 
cal time scales and a thermal relic of the Big Bang. Many 
scenarios for physics beyond the Standard Model (SM) 
include candidates for the DM particle, with a popular 
example being the neutralino in theories with weak scale 
supersymmetry (SUSY). 

Typically, in such theories the stability of the DM par- 
ticle is a consequence of a discrete Z2 parity symmetry 
motivated by the need to eliminate dangerous opera- 
tors allowed by SM gauge invariance. For example, in 
the Minimal Supersymmetric Standard Model (MSSM) 
renormalizable terms in the superpotential can violate 
baryon and lepton number and thus lead to rapid proton 
decay. R-parity conservation is often invoked to eliminate 
such operators, with the welcome side effect of rendering 
the lightest supersymmetric particle stable and therefore 
a DM candidate. 

It is remarkable that certain phenomenological con- 
straints such as proton decay may be overcome by invok- 
ing discrete symmetries, which in turn lead to stable DM 
candidates. However, our detailed knowledge of the fla- 
vor sector of the SM typically requires a highly non-trivial 
structure in the couplings of new particles to the quarks 
and leptons — one that is not easily achieved by imposing 
discrete symmetries. For example, R-parity does not pre- 
vent large flavor changing neutral currents (FCNCs) in 
the MSSM, and additional assumptions about the struc- 
ture of SUSY breaking must be made. In this regard 
the principle of Minimal Flavor Violation (MFV) is well 
motivated 0, One postulates that the global flavor 
symmetry of the SM is broken only by the Yukawa ma- 
trices in the presence of new particles. This provides a 
built-in protection mechanism against large flavor vio- 
lating processes. We shall demonstrate that, in analogy 
with the discrete symmetries discussed above, the MFV 
principle can lead to stable neutral particles and excellent 
DM candidates. 

In this paper we propose new theories of flavored DM. 



We add to the SM new matter multiplets with nontrivial 
transformation properties under the flavor group. We 
show that the MFV hypothesis implies that certain flavor 
representations contain one or more stable components, 
and we identify the general conditions for stability in 
these theories. Many variants of flavored DM models 
can then be constructed, each with a rich and predictive 
phenomenology. 

As an example, we study a simple effective theory of 
DM with a flavor SU{3)q triplet, gauge singlet scalar 
field, which provides a representative survey of various 
aspects of the cosmology and phenomenology in the fla- 
vored DM scenario. MFV predicts patterns in the mass 
splittings of the DM flavors governed by the SM Yukawa 
couplings. Flavor is intertwined with the physics of DM, 
as the Yukawa couplings and CKM matrix play a role in 
the thermal freezeout, direct and indirect signatures, and 
collider signals. In the scenario studied here, the heavy 
flavors in the DM multiplet decay to lighter flavors along 
with SM particles, leading to striking signatures at col- 
lider experiments. In other models, it is possible that the 
heavy flavors are stable on cosmological time scales, in 
which case several species comprise the cosmic DM. 

We note that several related studies exist in the liter- 
ature. As an alternative to R-parity, Ref. Q considers 
MFV as a protection mechanism in the MSSM, although 
in this scenario the LSP is not stable on cosmological 
times scales and is therefore not a DM candidate. Super- 
symmetric DM tied to the neutrino flavor structure has 
been considered in Q while Ref. Q investigated theories 
of flavored DM, though both of these proposals are out- 
side the MFV framework. Lastly, and in a more general 
context, models of flavor employing e.g. discrete symme- 
tries 0, Q or horizontal gauge symmetries P have also 
found their application in the DM problem. 

Flavor symmetries and dark matter. We be- 
gin by reviewing the MFV hypothesis, which provides 
a mechanism based on flavor symmetries by which new 
light particles may be consistent with precision flavor ob- 
servables. The quark sector of SM exhibits a large global 
flavor symmetry Gg = SU{3)q x 5'f/(3)„^ x 5[/(3)d„, 
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which is broken only by Yukawa interactions. When con- 
sidering theories of physics beyond the SM, the MFV 
principle dictates that new particles also respect Gq and 
the only sources of flavor breaking arise from insertions of 
the SM Yukawa matrices Yu and Yd [2, [3| • Formally, this 
is achieved by promoting the Yukawa matrices to spu- 
rion fields which transform nontrivially under Gq. The 
SM Yukawa interactions are written as 



Cy D QYddRH + QYuUbH^ + h.c. 



(1) 

(3,3,1), 



so that the Yukawa spurious transform as y„ 
Fd- (3,1,3) under Gq. 

A number of studies considering new MFV matter con- 
tent exist in the literature [10 15]. In all of these works, 
the Gq representations of the new particles have been 
chosen such that renormalizable tree-level couplings to 
the SM fermions are allowed. An immediate consequence 
of this is that such particles are unstable and decay to SM 
quarks or leptons. We now investigate the generality of 
this conclusion, and show in particular that MFV implies 
stability for particles in certain representations of Gq. 
With further assumptions about their electroweak quan- 
tum numbers, these stable particles become perfectly vi- 
able DM candidates. 

Let X be a matter multiplet that transforms non- 
trivially under flavor but is color-neutral, i.e. a singlet 
under SU{3)c. The most general operator that induces 
the decay of x then reads 

Cdocay ^ X Q ■ ■ -Q ■ ■ ■ Wfl^UR^ rfR^dfi^ (2) 

A B C D E F 

X Yu . . . Y^ . . . Yd . . . Y^ . . . Owoak, 

G ^~H^ I 

where A is the number of Q fields, B is the number of Q 
fields, etc. Cweak is a potential electroweak operator so 
that ^ is rendered invariant under SU{2) l x U{1)y. The 
flavor and color indices of the remaining fields are to be 
contracted by use of the SU{3) group invariant tensors 
6 and e so that ([2]) becomes a color and flavor singlet. 
This, however, can only be achieved if the triality ti of 
each SU{3)i tensor product {p,q)i with p factors of 3i 
and q factors of 3i, vanishes for the operator Odccay: 



U = {p~ q)i mods = 0, 



Q, ur, dji. (3) 



On the contrary, the decay operator © will be forbid- 
den if ti 7^ for at least one i. Denoting the irreducible 
representation of x under Gq as 

{nQ,mQ)Q X {nu,mu)un x {nd,md)dn, (4) 
where ng, mg, etc. can take values 0,1,2,..., the triality 



(n, m) 


SU{S)q X SU{S)u^ x SU{S)an 


Stable? 


(0,0) 


(1,1,1) 




(1,0) 


(3,1,1), (1,3,1), (1,1,3) 


Yes 


(0,1) 


(3,1,1), (1,3,1), (1,1,3) 


Yes 


(2,0) 


(6,1,1), (1,6,1), (1,1,6) 
(3,3,1), (3,1,3), (1,3,3) 


Yes 


(0,2) 


(6,1,1), (1,6,1), (1,1,6) 
(3,3,1), (3,1,3), (1,3,3) 


Yes 


(1,1) 


(8,1,1), (1,8,1), (1,1,8) 
(3,3,1), (3,1,3), (1,3,3) 
(3,3,1), (3,1,3), (1,3,3) 





TABLE I. Flavored DM candidates. Listed are the lowest- 
dimensional representations of Gq, organized according (n, m) 
where n = ng + Uu + rid, rn = mg + rriu + rud. We have also 
indicated the representations that are stable once MFV is 
imposed. Depending on their electroweak quantum numbers, 
these multiplets may contain viable DM candidates. 



conditions ([3]) become 

tc = {A-B + C-D+E-F)uYoA?. = Q, (5) 

tq = {nQ-mQ + A-B + G-H+I-J)iRod3^ 0, (6) 

tuj, = {nu-m.a + G-D-G+H)mod3^0, (7) 

tdn = {nd-md + E-F^I+ J) modi ^0. (8) 

Adding together Eqs. ([5][5]) and subtracting Eq. ([S]), we 
find that a necessary condition for Odccay to be allowed, 
and thus for x to be unstable is (n — TO)mod3 = 0, where 
n = uq + Uu + rid and m = mg -|- niu + rud- It follows 
that that Odccay is forbidden and x is stable if 



{n — m) mod 3 7^ 0. 



(9) 



Once © holds, x contains a stable component. In this re- 
gard, it is important to note that for operators with mul- 
tiple fields X, which may potentially mediate a loop 
induced decay, the stability condition above still holds. 
Table [T] lists the lowest-dimensional representations of Gq 
that are stable according to the condition 

In order to provide a viable theory of DM, any stable 
state in the fiavored multiplet x must further be electri- 
cally neutral. It thus remains to specify the electroweak 
quantum numbers of x- One possibility is that x is a SM 
gauge singlet [3] ■ Alternatively, as discussed in Ref . , 
X may be a n-plet of SU{2)l with hypercharge Y such 
that a component of x is neutral, Q = T^+Y = 0, with T3 
being the diagonal SU{2)l generator. According to this 
condition, a SU{2)l doublet has hypercharge Y — ±1/2, 
a SU{2)l triplet has hypercharge Y — 0, ±1, and so on. 

Note that the logic leading to the stability condition © 
made no assumption about renormalizability. Therefore, 
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the DM candidates are absolutely stable if we strictly en- 
force MFV, even at the nonrenormalizable level. We also 
made no assumptions about the Lorentz quantum num- 
bers of X) so one may consider scalars, fermions, vectors, 
and so on. Furthermore, if additional particle content is 
added besides the DM candidate our stability condition 
holds provided we demand that x contains the lightest 
mass eigenstate among the new particles. 

There exist variations in the precise implementation 
of MFV. For example, in the case of 'linear MFV, one 
imagines that each Yukawa insertion comes with a small 
coefficient, so that only the leading terms are important, 
whereas in 'general MFV all Yukawa insertions are con- 
sidered to be of similar size fisj . The number of Yukawa 
spurious in Eq. ([2]) was arbitrary, so that the stability 
condition ^ holds for general MFV. The case of linear 
MFV is more restrictive and may thus lead to additional 
DM candidates. 

Accounting for the different possible Lorentz, gauge, 
and flavor quantum numbers, we conclude that the MFV 
principle leads to a wide variety of DM candidates. MFV 
thus provides an attractive alternative to other well- 
known stabilization symmetries, such as the canonical 
Z2 parity, as well as higher Abelian or non-Abelian dis- 
crete or continuous symmetries [l9| . To be viable, these 
models must pass a number of constraints from cosmol- 
ogy and experiment. We now discuss these issues in an 
example model with a flavored DM candidate. 

Flavor SU{i)Q triplet dark matter. To illustrate 
in more detail the general phenomenological considera- 
tions in the flavored DM scenario, here we examine in 
some detail the physics of a simple flavored DM candi- 
date: a S'[/(3)q triplet, SM gauge singlet scalar fleld. 
We consider an effective theory with higher dimensional 
operators coupling the DM to quarks, and examine the 
constraints coming from cosmology, direct and indirect 
detection probes, and cosmology. We shall see that the 
theory has a rich phenomenology and leads to novel phe- 
nomena at hadron colliders. 

For concreteness, we consider a model with a scalar 
field with the following SU{i)c x SU{2)l x U{1)y x Gq 
quantum numbers: 



5~ (l,l,0)sM x (3,1,1)g, 
The Lagrangian is 



C 



V{S,,H)+C, 



(10) 



(11) 



where i = 1,2,3 is a flavor index, V{Si,H) is the scalar 
potential involving Si and the SM Higgs H, and Ccs 
contains higher dimensional operators involving two S 
fields and two quark fields. 

Spectrum. Let us first consider the scalar potential, 
which determines the mass spectrum of the flavored DM 
multiplet. The scalar S may have a bare mass term as 
well as contribution from electroweak symmetry breaking 



due to a Higgs portal coupling 16|. The relevant terms 
in the potential are 

V D mlSUa U, + b (Y^Y^)., +... )S, (12) 
+ 2A S:{a/1,, + h'{Y^Y^\, + . . . )S, H^H, 

where a, &, a', h' are dimensionless MFV parameters. The 
ellipsis indicate further MFV spurion insertions involving 
Yd which generally lead to mass fine-splittings; for sim- 
plicity we neglect these terms. The Higgs field obtains 
a vacuum expectation value {H) = u/-\/2, breaking the 
electroweak gauge symmetry and giving a contribution 
to the S mass-squared matrix: 

c D -s* [m^i., + m|(y„r„t),,] S,, (13) 

where we have defined 

2 2 I \ 2 / 
rriA = TOoa -f Af a , 



lib- 



■2U 



(14) 

Without loss of generality, we freeze the Yukawa spu- 
rious to the background values Yd — \di Yu = V^X^ 
where (A„)y = yu,A], {^d)tj yd.t^tj are diagonal in 
the physical Yukawa couplings yu^i and ydX, quark mass 
eigenstates are then obtained via ul — J> V^'ul. The mass 
matrix in Eq. ([T5)) is diagonalized by 5 — >■ S: 



-S* [mil 



5, 



(15) 



Up to fine splittings induced by up and charm Yukawa 



couplings, we find mf 



and 



. There are two possibilities for the spectrum: 1) 
a "normal" spectrum with two lighter, nearly degenerate 
states 51,2 and one heavier state S'3, or 2) an "inverted" 
spectrum with two heavier states 5*1.2 and one lighter 
state S'3. 

Of phenomenological importance are the trilinear cou- 
plings of scalars to the Higgs h jl6|, which in the mass 
eigenbasis read: 



C^-2\vhS*{a'U,+b\\l),,)S, 
= -~2XivhS* Si, 



(16) 



where we have defined A^ = A(a' + h'y^ .j) These cou- 
plings provide annihilation channels into SM particles, 
and are furthermore constrained by direct and indirect 
DM searches. For light DM particles in the 10 - 100 
GeV range, the recent XENONIOO nuU resuhs ^ im- 
ply a limit of Xi < 0.015 (m,/10 GeV)(m/i/120 GeV)^ 
for the DM particle. In addition, the Higgs portal inter- 
actions above lead to nonstandard decay channels of the 
SM Higgs, h — >• S*Si, provided that the scalars Si are 
light enough. The partial width for h — > S* Si is given by 



1 



Ami 



1/2 



(17) 



For values of Xi on the order of the 5-quark Yukawa cou- 
pling yb, these decay modes become competitive with the 
SM h ^ hh mode for a light Higgs. 
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Effective theory. By construction, the scalars do not 
couple to quarks at the renormahzable level since such 
operators lead to the decay of the would-be DM can- 
didate. However, couplings between pairs of DM multi- 
plets and quarks may occur in an effective theory through 
higher dimensional operators. At the dimension six level, 
we can write the following effective operators coupling 
two scalar multiplets to a quark and anti-quark: 

ol.ki = {Q^rQ,){SlKsl). (18) 

OIm = {umi^un,){S*Xs,), (19) 

0%M = idmrdR^XS-lttSi), (20) 
Ot,k, = {Q^UR,)iS;Si)H^ + h.c, (21) 
Ol^,, = iQ,dR,){S:Si)H + h.c., (22) 

where i, j, k, I are flavor indices. From these operators we 
can construct the effective Lagrangian: 

1 ^ 

•^e// = -Tj ^ (^ijki^ijki- (23) 
1=1 

The coefficients cj^f.^ in general contain all possible flavor 
contractions consistent with MFV. For example, for the 
operator we have 

(r„r„t)« + ci(y„y„t) 

+ cl*lu{Y^Y^)k, +■.., (24) 

where c] are Wilson coefficients for a given flavor con- 
traction, and where we display the leading terms in the 
MFV expansion up to one insertion of YuY^. 

Such higher-dimension operators can have a signiflcant 
impact on the cosmology by providing the DM with effi- 
cient annihilation channels to SM quarks, as long as the 
effective scale A/ ^/c^ suppressing the operators is not too 
much larger than the weak scale. The same operators can 
potentially be probed by direct and indirect detection ex- 
periments, colliders, and precision flavor studies. 

Indeed, certain operators and flavor structures are 
more constrained than others. To illustrate this fact, 
consider for concreteness the "inverted" spectrum from 
Eq. (|15|) in which the DM particle is S3 . For the operator 
O^, the flavor structures associated with the coefficients 
cj and C4 in Eq. lead to the vector coupling of the 
DM particle S3 to valence quarks in the nucleon, and are 
therefore strongly constrained by direct detection exper- 
iments. On the other hand, the flavor structures associ- 
ated with C2, C3, and C5 imply that the vector couplings 
of the DM state 5*3 are dominantly to the third genera- 
tion quarks, while the couplings to valence quarks either 
vanish or are Yukawa and/or CKM suppressed. There- 
fore, these coefficients are far less constrained by direct 
DM searches. 

Continuing on with this example, there is the possibil- 
ity of new FCNCs induced by the operator once the 



DM multiplet is integrated out. Whereas the coefficients 
cj and c\ in Eq. p4)) lead to strictly flavor-diagonal cou- 
plings in the quarks, the flavor structures associated with 
C2, C3, and C5 indeed do induce FCNCs in the down-quark 
sector due to the presence of CKM mixing. 

A number of other phenomenological considerations 
must be taken into account, each of which potentially 
constrain the Wflson coeffcients cf^^.; for the operators in 
Eqs. (I19ll22p . Depending on the UV completion some of 
these coefficients may be sizable while others may be sup- 
pressed or vanish. From the effective theory perspective, 
one can explore the consequences of each operator on the 
cosmology and phenomenology in a model-independent 
way. Understanding which operators are constrained and 
which are allowed may then point towards the underly- 
ing dynamics at scales greater than A. A comprehensive 
analysis of the effective operators (I18II22P and their vari- 
ous flavor structures will be left for future work. Instead, 
we shall specialize to one particular operator and examine 
in detail the associated cosmology and phenomenological 
constraints and prospects. 

For the remainder of the paper let us consider operator 
in Eq. ((22|) with the following flavor structure: 

Caff = ■^[Q^S,][S*{Yd)JkdRk]H + h.c, (25) 

where the flavor indices are contracted within the brack- 
ets, while the Lorentz, color, and SU{2)l, indices are 
contracted outside. The coefficient c is taken here to be 
a phase, so that |c| = 1. As we will demonstrate, this 
operator allows for a viable cosmology while being con- 
sistent with a variety of constraints. 

After electroweak symmetry breaking and diagonaliza- 
tion of the quark and scalar masses, we are left with the 
following effective Lagrangian: 

C,ff ■^^[dLrV^^S,][SliVXd)kedRi] + h.c. . (26) 

Focusing on the "inverted" spectrum below, nij < mi ~ 
1712, S3 is the stable DM candidate, while Si and 5*2 are 
unstable. 

Relic abundance. First we consider the thermal relic 
abundance of 5*3 . If the DM mass 777.3 is larger than the 
bottom quark mass rUb, the dominant annihilation mode 
is ^a^g — ?> bb, coming from the following term contained 
in Eq. 

C,s D -^mhlVtbl^S^SsbLbR + h.c. , (27) 

where |Vtf,| ~ 1 is the top-bottom CKM matrix element. 
The above term in the Lagrangian leads to the following 
thermally averaged annihilation cross section: 

(--)33^S. = 4;^-?l^*.l'(l-^)'^' (28) 
x|[Re(c)p (^l-!!|)+[Im(c)]^}. 
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In the limit ma ^ uib, the annihilation cross section 
approximately scales as 



(av) 



33-^bb 



^ , /200GeV\^ 



(29) 



suggesting that an effective scale of A « 200 GeV is re- 
quired in order to match with the nominal 3a DM abun- 
dance il^j^jh^ — 0.1109±0.017 as inferred from the seven 
year data sample of the WMAP satellite experiment f2l| . 
This can also be seen in Fig. [T] where the gray shaded 
band indicates the region in which 6*3 -I- ^3 attain the 
right relic abundance. Below the threshold of bb an- 
nihilation, 7713 < 4.2 GeV, the most efficient channel is 
SIS3 ^ sb,bs with (av) cx \Vts\^\Vtb\^mbm,/A\ 

CMB constraints. Residual annihilation of S'3-pairs 
has its strongest impact at cosmic redshifts z ^ 1000 
at the time of recombination of electrons with hydrogen 
(and helium.) The electromagnetic cascades formed in 
the annihilation delay recombination thus affecting the 
temperature and polarization anisotropics of the cosmic 
microwave background radiation (CMB) [S^,!!^!- Since 
the rate of injected energy scales as {crv)n'jjj^jmDM oc 
(av) I m DM where hjjm is the DM number density, the 
CMB limit on the annihilation cross section can conve- 
niently, and in a rather model-independent fashion, be 
expressed in the form 



(1 - fi^){<7v)/m3 < r. 



(30) 



Here, fn is the fraction of energy carried away by neutri- 
nos and we employ r = 0.191 GeV/cm'^s"^ as obtained 
in Using the PYTHIA [23| model for e+e" annihila- 
tions to liadrons we have verified that always remains 
below 10%; as a representative value we choose = 0.08. 
The resulting lower limit on the effective scale A is shown 
by the dashed line in Fig. [T] and is strongest for small val- 
ues of 7713 , excluding the otherwise cosmologically viable 
region < 5 GeV. 

Direct detection. The effective operator p6|) also in- 
duces the elastic scattering of 5*3 on nuclei. The following 
terms lead to spin- independent scattering: 



Re(c) 



3 

E 

1=1 



\V3,\'S^SM 



(31) 



Though the couplings have the usual dependence on the 
quark masses, as for example expected in Higgs-mediated 
DM scattering on nuclei, the couplings to light quarks s 
and d are additionally CKM suppressed. This implies 
that only the scattering of 5*3 on the 6-quark content of 
the nucleons is relevant in constraining the scale A as well 
as the phase c of the higher-dimensional operator ([26| . 
The nucleon matrix element under question reads 



fn,b = {n\mbbb\n) = m, 



2 

'27 



„(n) 

■Itg 



0.04, 



(32) 



lO'^ 



> 

o 



10^ 



10^ 




ms [GeV] 



FIG. 1. Constraints; Shown are lower limits on the eflFective 
scale A as a function of ms for Re(c) = 1. The dashed line is 
a constraint on the annihilation cross section from CMB tem- 
perature and polarization anisotropics. The solid (red) line 
is the exclusion bound from the recent XENON 100 direct de- 
tection experiment with 100 live-days of exposure; the dotted 
(blue) line shows the sensitivity of a future ton-scale liquid 
xenon experiment. In the lower right shaded region, the ef- 
fective field theory description breaks down. Within the gray 
shaded band, the relic density of 5*3 matches the observed DM 
density, ^s^+s^h'^ 



where we have used that the fact that the 6-quark matrix 
element can be expressed in terms of the matrix element 
of the gluon scalar density fisj, /^^ — 1 — X^u d s fr'q — 
0.57 [26,]. 

This leads to the effective spin-independent elastic 
DM-nucleon cross section. 



mcr\vtb?fl,^i: 



(33) 



43 2r , ,,2 /10GeV\ VsOOGeVX"* 
.3xlO-cm^[Re(c)f (^-— j (^^^ j . 

The upper limit on (t„ from direct detection experi- 
ment can be translated into a lower limit on A. The 
most constraining data in this respect comes from the 
recent 100 live-day result of the XENONIOO collabora- 
tion [20] . We obtain the constraint by using Yellin's max- 
imum gap method [27| which accounts for the 3 observed 
events within the 48 kg fiducial detector volume and em- 
ploy the results of [28| to account for detector resolu- 
tion, efficiency, and acceptance. The resulting limit in 
the {rnnM, fn) plane is in very good agreement with the 
one of [20|. For the sake of exploring future sensitivity 
of direct detection experiments we also simulate a ton- 
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scale liquid xenon detector assuming a raw exposure of 
1 ton X year using the detector model for XENONIO of 
Ref. H. 

Choosing Re(c) = 1, the XENONIOO lower limit on 
A as a function of is shown by the solid (red) line 
in Fig. [H As can be seen, mass values in the range 
10 GeV < ma < 60 GeV are excluded. The dotted (blue) 
line illustrates that a future ton-scale direct detection ex- 
periment may well be sensitive to the entire electroweak- 
scale parameter region of m^, excluding values up to 
ma < 400 GeV. It is important to note, however, that 
since the overall phase c of the higher dimensional op- 
erator is a free parameter, one can avoid completely the 
direct detection constraints by choosing Re(c) ^ 1 while 
still obtaining the correct relic abundance. Finally, the 
the shaded region in the lower right corner of Fig. [1] — 
subject to the condition A < m3/27r — indicates the part 
of the parameter space that does not admit a perturba- 
tive UV-completion [29{ . 

Flavor physics. The new scalars can lead to new FC- 
NCs at the one loop level. Here we will focus on the 
contribution to K^~K^ mixing. The relevant interac- 
tions from the operator ()26p that can in principle induce 



mixing are 



(34) 



Integrating out the scalars one arrives at the following 
Lagrangian: 



Cf{sRdL){sRdL) + h.C. 



where the Wilson coefficient is given by 



CI 



327r2A4 



'■F[ ^ 
mi 



(35) 



(36) 



withF(x) = (5: + l)(.x-l)-ilogx-2. The function F(i) 
takes 0(1) values for typical values of the masses, and 
vanishes for x = 1. To arrive at Eq. ([36]) we have made 
use of the unitarity of the CKM matrix and assumed a 
mass splitting predicted by MFV With A = 0{y) as 
dictated by the relic abundance calculation, we observe 
that the operator ((26|) leads to a contribution in Eq. ((36)) 
which is suppressed compared to the SM value by a factor 
rn^/rriY/ ~ 10~^. Therefore, we are comfortably below 
the experimental limits. Similarly, for Bq — Bq mixing, 
we find operators that are suppressed compared to the 
SM value by a factor rn^/rn^ k, 10^'^. Furthermore, we 
find that no contribution to the inclusive decay b ^ sj 
is induced. 

Decays of heavy dark flavors. The heavier scalars Si 
and 15*2 are unstable and decay via the higher dimen- 
sional operator (pS)) . The leading operators which are 
unsuppressed in powers of light quark masses ms.d and 



off-diagonal CKM elements are 

CcS 3 j^mbVtbV*^S^S2SLbR 



j^mbVtbV:aS;SidLbR + h.c. . (37) 



These operators mediate the decays 

5*2 S3sb, 
Si S^db. 



(38) 
(39) 



In the limit mi 2 3> mb, m^, the three-body partial decay 
width for Si — >■ S^qb is approximately given by 



i—¥3qb 



\cVtbV*\'mim: 
5127r3A4 

lO^^MeV X 



(40) 



25 GeV 



200 GeV y 
A j 



The scalars therefore decay promptly. 

LHC prospects. A novel prediction of our fiavored DM 
scenario is the presence of heavy dark flavors. Given a 
production mechanism, such heavy flavors can lead to in- 
teresting signatures at the Large Hadron Collider (LHC). 
The effective operator (j26p can in principle lead to pair 
production via the parton-level processes qiqj — SiS*. 
In practice, however, the production cross section for 
these processes is negligible due to the Yukawa sup- 
pressed couplings in Eq. ([^ (this is also why the oper- 
ator faces no bounds from monojet searches at the Teva- 
tron 

Here we instead focus on production of Si^2 via the 
Higgs portal interaction in Eq. (jl3l) . If the new scalar 
states are lighter than m/i/2 and have a trilinear coupling 
to the Higgs in Eq. (|T5)) larger than the bottom Yukawa 
interaction yb, then the Higgs decays predominantly to 
intermediately heavy scalars 5*1^2 with a partial width 
given by Eq. ([T7)) . These scalars subsequently decay via 
the three-body processes (|38l39p to a bottom quark, a 
light quark {d or s), and the DM candidate 5*3, 



h S*iSi bbddS;S3, 
h S^S2 bbssS^S 



(41) 
(42) 



We now investigate the potential of the LHC {y/s — 14 
TeV) to discover the processes (j41l42p for the benchmark 
scenario given in Table [Hi For this choice of parame- 
ters we find that the Higgs decays almost exclusively via 
(|41|42|) . The most promising channel to disentangle the 
signal from the large SM backgrounds is expectedly the 
HZ production channel, where the Z boson subsequently 
decays to electrons or muons. One immediate reason for 
this is that the resulting hard leptons provide a trigger- 
able signal. Furthermore, if the Z boson is boosted and 
the Higgs recoils against the leptonically decaying Z, the 
leptons yield a good handle to calibrate the jet energy 
and to measure the amount of in the event. Similar 
configurations have been used in previous analyses for 
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mh 


mi 


1712 


ma 


Ai 


A2 


As 


m| 


A 


a 


a 


6 


b' 


120 


25 


25 


10 


0.15 


0.15 


0.01 


-86^ 


0.35 


1.14 


0.429 


-1.10 


-0.409 



TABLE II. Benchmark spectrum and couplings for LHC study. The input parameters are mi, m2, ma, Ai, A2, A3 
This determines the Lagrangian parameters a, a , b, b' . All mass parameters are in GeV. 



m| and A. 



a boosted SM Higgs boson which decays into a bottom 
quark pair [s^ . 

We simulate the signal using MadEvent jsij and the 
backgrounds using Herwig-|--|- [32]. The NLO production 
cross sections for the backgrounds and the signal was cal- 
culated using MCFM |3j| . To select a lepton we demand 
it to be central and sufficiently hard. 



<2.5, 



PT.i > 20 GeV. 



(43) 



For the hardest lepton in the event we further require 

PT,h > 80 GeV. (44) 

The niuons and electrons must be isolated, i.e., the ha- 
dronic transverse energy in a cone of R = 0.3 around 
the lepton has to be i^Thadro„ic < 0.1 Et.i. We accept 
events with exactly two isolated leptons. The Z bosons 
are reconstructed by combining two oppositely charged 
isolated electrons or muons, requiring 



mz - 10 GeV < mu < mz + 10 GeV 



(45) 



and pt,u > 150 GeV. 

We group the remaining visible final state particles into 
"cells" of size Aij x Acjj = 0.1 x 0.1. All particles in a cell 
are combined, and the total three-momentum is rescaled 
such that each cell has zero invariant mass. Cells with 
energy < 0.5 GeV are discarded, while the remaining 
ones are clustered into jets. We recombine the jets using 
Fastjet 3JI with the anti-kT jet algorithm 35|, i? = 0.7, 
and 

PTj > 30 GeV. (46) 

We reject events where a jet with prj > 50 GeV has a 
smaller angular separation to the Z boson than Ai?^ ^ < 
1.5. At least one jet must be present in the event, and 
the hardest jet in the event must be 6-tagged. For the b- 
tag we assume an efficiency of 60% and a fake rate of 2%. 
The 5*3 states escape the detector unobserved. Therefore 
we accept an event if 



> 50 GeV, 



(47) 



and calculate from all visible objects with \y\ < 5.0 . 
The resulting vector must have a large azimuthal 
separation from the reconstructed Z boson. 



FDM 
tt 

VV tTTTT), 




50 100 150 200 

mT(S(i,2,,S(i.2)) [GeV] 



300 



FIG. 2. Reconstruction of mr according to Eq. (|49|) for a 
Higgs mass of 120 GeV. In addition to the cuts outlined in 
Eqs. (|43p - (|48p we require at least two jets to accept an event. 



With these cuts we obtain S/B ~ 1 and S/ \fB ~ 5 after 
15 fb~^; see Table Hill for details. Therefore an excess in 
events with J^t + bb + can be observed relatively 
early at the LHC. However, due to the sizable fraction 
of missing transverse energy, fJr, the reconstruction of 
the mass of the Higgs boson is challenging. Because the 
bottom quark and the light quark from the three-body- 
decay of S'1,2 are highly collimated they will be observed 
as one jet. Thus, after requiring at least two jets in the 
event, in the rest frame of the Higgs one finds m^^Sj ~ 
rujj . Following [s^ , we use 

tot(5'*S'j) = 



o^-(PT,,,■+I^T)^ (49) 



to reconstruct the Higgs mass as shown in Fig. [5] In 
Eq. (|3ni)i the transverse energies are given by = 

-I- m^,- and Etjj — 



rrijj is the in- 



A(j)p^,z > 2.0. 



(48) 



variant mass of the two hardest jets in the event which 
have a large angular separation to the Z boson, Rj z > 
2.0. For a rough mass reconstruction of the Higgs boson 
in this channel we expect 0(100) fb""'^, but other Higgs 
production channels could be added. If mi^2 ^ W3 -I- mi, 
the Higgs decays to jets and missing transverse energy. 
This signature displays strong similarities with the so- 
called buried Higgs models [37|. In either case, jet_sub- 
structure can help to reconstruct the Higgs boson ]3B\ . 

Discussion. In this paper we have demonstrated 
that the MFV hypothesis, motivated by the need to sup- 
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FDM 


tt 






WW 


Tij > i, ni — I and Pt,Ii > oU OeV 


1 O ^7 


o9Uo. / 


202. 6 


ioo.o 


242. z 


> 50 GeV 


7.8 


5744.1 


20.6 


20.4 


118.8 


Z reconst. and pt,z > 150 GeV, no ARj^g^z < 1-5 


4.3 


9.9 


5.8 


3.8 


0.7 


A<?!>9j,,z > 2.0 


4.2 


4.6 


5.2 


3.3 


0.03 


b-tag 


2.2 


2.2 


0.2 


0.1 


0.01 



TABLE III. Cross sections for different steps in the analysis for the signal from the flavored dark matter model (FDM) and the 
dominant SM backgrounds. Systematic uncertainties in the measurement of can elevate Z+jets to a relevant background 
but their simulation is beyond the scope of this work. All cross sections are normalized to NLO QCD precision, calculated using 
MCFM, and are given in femtobarn. The background events were generated using Herwig++ and the signal was generated 
using MadEvent and showered using Herwig++. 



press dangerous flavor changing processes in new physics 
scenarios, also provides a novel organizing principle for 
the physics of DM. Indeed, MFV leads to the stability of 
the lightest neutral particle for a large number of flavor 
representations. We have determined which representa- 
tions have stable particles. This opens up the possibility 
for a large class of new flavored DM candidates. 

Because of the connection with flavor physics, the phe- 
nomenology of these models is very rich. In addition 
to providing stability, MFV is a predictive framework 
for DM. The mass spectrum of the flavored DM multi- 
plet is governed by the hierarchical nature of the quark 
Yukawa couplings. Furthermore, the pattern of couplings 
between the DM and the SM is prescribed by the fla- 
vor symmetries of Gq and CKM breaking. Our example 
study of the effective theory of a flavor SU{3)q triplet, 
gauge singlet scalar DM illustrates many of the general 
aspects of the cosmology and phenomenology, and leads 
to a viable DM scenario with testable predictions at fu- 
ture direct detection experiments and present day collid- 
ers. In particular, the signatures of the heavy unstable 
dark flavors are quite novel and can be searched for at 
the LHC. However, it is fair to say that we have only 
scratched the surface of the phenomenological implica- 
tions of the flavored DM scenario. 

In this paper we have focused on the flavor symmetries 
of the quark sector, primarily because of our still incom- 
plete knowledge about the SM neutrino sector. How- 
ever, with some assumptions about how neutrinos ac- 
quire mass the MFV hypothesis can be formulated for 
leptons [39], and it would therefore be interesting to in- 
vestigate the possibility of leptonically flavored DM. In 
this case, it is conceivable that heavy dark flavors may 
leave their imprint via decays to DM with associated lep- 
tons. 

It would be interesting investigate more generally the 
predictions for other flavored DM multiplets, including 
larger flavor representations and electroweak multiplets. 
For the latter, typically because of the efficient annihi- 
lation to SM gauge bosons, the DM particle must be 
quite heavy to be in accordance with the observed cosmic 
abundance 17|. Multiple flavors of stable electroweak 



DM may instead be comparatively light since each flavor 



then only constitutes a portion of the DM density. 

More comprehensive studies of flavored DM effective 
theories is certainly warranted. Much attention as of late 
has been given to effective fleld theories of DM, notably 
as a means to investigate the potential of hadron collid- 
ers to probe DM-quark interactions in a manner comple- 
mentary to direct detection experiments [29]. In these 
studies, the couplings of DM to quarks were taken to 
be flavor universal for vector couplings and proportional 
to mass for scalar couplings to avoid FCNC constraints. 
Of course, a less restrictive possibility is to abide by the 
MFV principle, which allows for more general couplings 
between DM and quarks, including flavor violating cou- 
plings governed by the CKM elements. 

Ultimately, since flavored multiplets cannot couple to 
quarks at the renormalizable level (else they would de- 
cay) , it is interesting to construct UV completions to such 
effective theories of flavored DM. One may postulate ad- 
ditional particles that connect the DM multiplet to SM 
quarks, which themselves must also obey MFV. Indeed 
it is certainly interesting to speculate about a new sector 
of MFV matter, given that MFV has proven useful in 
addressing recent Tevatron anomalies, such as the like- 
sign dimuon asymmetry in B-decays and the top 
quark forward-backward asymmetry 42. |43|. Perhaps, 
the DM particle is contained in such an MFV sector. 

The evidence for DM is one of many reasons to expect 
new physics, and yet the SM remains impeccable when 
confronted with experiment. This motivates new symme- 
try principles designed to forbid dangerous interactions 
of new particles, which may also have the dual purpose of 
stabilizing DM. The canonical example of such a symme- 
try is a discrete Z2 parity. In this paper, we have shown 
that flavor symmetries together with the MFV principle 
provide another compelling possibility, thus pointing the 
way towards new avenues for physics beyond the SM and 
DM. 
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